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We study two models of a system of queues served in cyclic order by 
a single server. In each model, the \th queue is characterized by general 
service time distribution junction H;(-) and Poisson input with param- 
eter X s . 

In the exhaustive service model, the server continues to serve a particular 
queue until the server becomes idle and there are no units waiting in that 
queue; at this time the server advances to and immediately starts service 
on the next nonempty queue in the cyclic order. 

The gating model differs from the exhaustive service model in that when 
the server advances to a nonempty queue, a gate closes behind the waiting 
units. Only those units waiting in front of the gate are served during this 
cycle, with the service of subsequent arrivals deferred to the next cycle. 

We find expressions for the mean number of units in a queue at the 
instant it starts service, the mean cycle time, and the Laplace- Stieltjes 
transform of the cycle time distribution function. 

I. INTRODUCTION 

We consider a system of queues served in cyclic order by a single 
server. The tth queue is characterized by general service time dis- 
tribution function #((•) and Poisson input with parameter A*. 

We study two variations of this model. In the first, called the 
exhaustive service model, the process begins with the arrival of a 
unit at some queue, say A, when the system is otherwise empty. The 
server begins on this unit immediately, and continues to serve queue 
A until for the first time the server becomes idle and there are no 
units waiting in queue A. The server then looks at the next queue 
in the cyclic order, queue A + 1, and serves those units, if any, that 
have accumulated during the serving period of queue A. The server 
continues to serve queue A + 1 until for the first time the server 
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becomes idle and there are no units waiting in queue A + 1. The 
process continues in this manner, with the queues being served in 
cyclic order, until for the first time the system becomes completely 
empty. The process is then re-initiated by the arrival of the next 
unit. No time is required to switch from one queue to the next. 

The second variation, called the gating model, differs from the first 
in the following way: When the server moves to a queue with at 
least one waiting unit, the server accepts only those units that were 
waiting when the server arrived, deferring service of all subsequent 
arriving units until the next cycle. That is, in the gating model, at the 
instant the server advances to a nonempty queue a gate closes behind 
the waiting units, and only those units waiting in front of the gate are 
served during that cycle. 

The exhaustive service model is analyzed in detail. We obtain the 
generating function for the joint probability distribution of the num- 
ber of units in each queue at an instant at which the server finishes 
serving any one of the queues. We then obtain expressions for the 
mean number of units in a queue at the instant it starts service, the 
mean cycle time, and, in a form suitable for numerical computation, 
the Laplace-Stieltjes transform of the cycle time distribution function. 

Finally, we note that the equations describing the gating model 
differ only trivially from those of the exhaustive service model, and 
that the same method of solution applies to each. 

Systems in which a single server is shared among several queues 
are common. For example, in the No. 1 Electronic Switching System 
the central control spends much of its time polling various hoppers 
and performing work requests that it finds in these hoppers. Similarly, 
in a time-shared computer system the users have access through tele- 
typewriters to a central computer which is shared among them. The 
cyclic queueing models studied here are of a type which may be useful 
in the analyses of these and similar problems. 

The exhaustive service model for the special case of two queues 
has been studied by L. Takacs, 1 B. Avi-Itzhak, W. L. Maxwell, and 
L. W. Miller, 23 and M. F. Neuts and M. Yadin.* Avi-Itzhak et al 
used an argument based on the properties of mean values, to obtain 
an expression for the mean waiting time suffered by a unit in either 
queue. Takacs, by a more direct argument, utilizing the Markov 
chain imbedded at the epochs of service completion, obtained the 
corresponding Laplace-Stieltjes transform and formulas for the wait- 
ing time moments assuming service in order of arrival. Neuts and 
Yadin obtained waiting time results for the transient case. 
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The more general two-queue model in which the time required to 
switch from one queue to the next has some arbitrary distribution 
function is also studied in Ref. 3, and has been investigated in addi- 
tion by M. Eisenberg 5 and J. S. Sykes. 6 M. A. Leibowitz 7 - 8 has stud- 
ied a multiqueue model similar to the gating model studied here. A 
nonprobabilistic approach to cyclic queueing problems has been used 
by J. B. Kruskal. 9 

In the present paper we use the imbedded Markov chain approach 
but, as with Neuts and Yadin, our chain is imbedded at the instants 
at which the server completes serving a queue, rather than at the set 
of all instants of service completion used by Takacs. Whereas Takacs 
and Neuts and Yadin obtained waiting time results, our analysis yields 
cycle time results. The mathematical analyses characterizing the three 
approaches share some common ground, although the differences, espe- 
cially those arising from our consideration of an arbitrary number of 
queues, are significant. 

Also, in a recent nontechnical article on queues by Leibowitz, 10 the 
present problem is offered as a prime example of an important, dif- 
ficult, unsolved queueing problem. 



II. PRELIMINARIES 



In the analysis of the exhaustive service model, we take the number 
of queues to be N + 1 ^ 2. Units arrive at the ith queue according to 
the Poisson process with rate X, ; that is, the probability Q,-(fc; t) that 
k units arrive at the ith queue in an interval of length t is 

Qi (k;t) = ^exp(-M) (ft « 0,1,2, ■•• ;i = 0,l, ••• ,N). 

The length of time required to serve a unit from queue i has distribution 
function //,-(•) with mean h { (i = 0, 1, • • • , N). 

In the analysis of the exhaustive service model, we shall use the 
concept of busy period, discussed at length by Takacs 11 . For the ordinary 
single-server queue, the busy period is defined as the length of time 
from the instant a unit enters a previously empty system until the 
next instant at which the system is completely empty. Both the distri- 
bution function of the busy period and its Laplace-Stieltjes transform 
are known explicitly for the M/G/l queue. In particular, the M/G/l 
queue with arrival rate X and mean service time h has a busy period 
with mean b = h/(l - \h) if \h < 1 and b = » if \h ^ 1. 

Consider now the M/G/l queue with j waiting units; define the 
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j-busy period as the length of time from the instant at which service 
starts on the first of the j units until the next instant at which the 
system is completely empty. (When j = 1, the j-busy period and the 
busy period are identical.) Each of the j units, which together generate 
a ;-busy period, individually generates a 1-busy period. Thus (as 
Takacs shows) the distribution function of the j'-busy period is the 
j'-fold convolution with itself of the distribution function of the 1-busy 
period. 

Denote by B { (-) the distribution function of a 1-busy period for 
queue i, by /3,-(-) its Laplace-Stieltjes transform, and by 6, = hj{\ — 
\ihi) its mean. Let £?'(•) be the ;'-fold convolution of !?,•(•) with itself, 
B*}(0 = Bi(-). Then a i-busy period for the t'th queue has distribution 
function Bf'(-) and Laplace-Stieltjes transform (£«(•))'■ 

III. FORMULATION OF IMBEDDED MARKOV CHAIN STATE EQUATIONS FOR THE 
EXHAUSTIVE SERVICE MODEL 

There are N + 1 ^ 2 queues. Suppose that the system is idle and a 
unit arrives at some queue at epoch t . The server immediately com- 
mences service at that queue, and continues to serve units at that 
queue until the first instant r l at which that queue becomes empty. 
If the system is not empty at r x , the server advances to the next queue 
in the cyclic order. The server immediately commences work at this 
queue until the instant t 2 at which this queue becomes empty (where 
T2 = Tl if the server finds the queue empty), and continues on in this 
manner until for the first time, t„ say, the server finishes serving a 
queue and there are no units waiting anywhere in the system. The 
process terminates at r n and is reinitiated by the next arrival. 

Thus the process generates a set of points t , tj , • • • , t„ , where t 
is the arrival instant of a unit at some queue in the previously empty 
system, and r n is the first instant at which the system becomes com- 
pletely empty again. The next arrival, at epoch r » say, reinitiates the 
process, and a new set of points, t - , t v , • • • , t„. is generated. We call 
the points t x , • • • , t„ (and t v , • • • , r n .) switch points. 

Note that r is not a switch point, whereas t„ is a switch point. 
Successive switch points may occur simultaneously in time, but are 
nevertheless considered distinct. Thus, with each switch point is as- 
sociated a queue, namely, that queue at which the server has just 
completed its visit. 

When the server finishes serving a queue and finds the system com- 
pletely empty, a switch point associated with that queue is recorded. 
The next switch point is recorded when the server leaves the queue at 
which the process is reinitiated, and is associated with that queue. 
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Let (i; n Y , • • • , n N ) denote the state of the system at an arbitrary- 
switch point, where i is the index of the associated queue, and n k is the 
number of units waiting in queue i + k(k = 1, • • • , N). [For simplicity, 
no special notation will be used to denote arithmetic mod (N + 1).] 
Let the state (i; n t , • • • , n N ) have probability Pi(n x , • • • ,n N ); that is, 
P,(?ii , • • • , n N ) is the joint probability that at a switch point, the 
server has just completed a visit to queue i (i = 0, 1, • • • , N) and ?ij 
units are waiting in queue i + 1, n 2 units in queue 2 + 2, • • • , and n N 
units in queue i + N. 

The state (i; nj , • • • , n N ) can occur through the following exhaustive 
and mutually exclusive contingencies: 

(i) The server leaves queue i — 1 and finds j ^ 1 units waiting 
for service in queue i, where it thus spends a length of time 
equal to a j'-busy period. 

(it) The server leaves queue i — 1 and finds j = units waiting 
for service in queue i, but at least one unit waiting for service 
somewhere else in the system, so that the server then 
"passes through" queue i in zero time. [That is, the state 
(i; n l , ••• %_, , 0) necessarily follows the state (i — 1; 0, 
fti >• • •> ftjv-i) where at least one of the n k 9 £ 0(k = l,- • •, iV— 1).] 

(in) The server leaves some queue and finds no units waiting any- 
where in the system. With probability X.-/X (X = X + • • • + X^) 
the next arrival (which reinitiates the process) occurs at queue i, 
where the server then spends a 1-busy period. 

These considerations lead directly to the imbedded (at the switch 
points) Markov chain probability state equations: 

P,(rii , • • • , n jV ) 

= E Z • • • £' Pi-iU, fci . • ■ ■ . **-0 

,=) k,-0 t»-;=0 

• r II Qi+ m (n m - k m ; t)Q^ N (n N ; t) dBf(t) 

Jo m-1 

+ P,_,(0, n, , • • • , n iV -!) (l - s( Z n„jj 8(n N ) 
+ vt P 4 (0, ■■ ,0) T fl Q i+m (n m ; /) dB,(t) 



il if 1 = • n i *T 

8(x) = 1 ; i = 0, I, ■ ■ ■ , N 

.0 if x ^ 



(1) 
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[Throughout the analysis, arithmetic mod (N + 1) in subscripts 
will not be specially denoted.] Assuming it exists, the distribution 
{Pi(rii, • ■ • , n N )} is uniquely determined by (1) and the normaliza- 
tion equation 



If oo 

— ^ 

-0 



E E ••• £P«(n. , ••• ,n N ) = 1. (2) 



(Intuitively, one would expect a unique stationary distribution to 
exist when E*--o A A- < !•) 

IV. FUNCTIONAL EQUATIONS FOR GENERATING FUNCTIONS 

We define the probability generating functions gi(x±, • ■ • , ##) : 

g t (x t ,*■•,»*■)= 2 ■■• £ ^<(ni , • • • , w^s?' • • • x^ w 

(»- 0, 1, ... ,JV). (3) 
Substitution of (1) into (3) yields, after some rearrangement, 

gt(xi , • • • , x„) 

- S £ ■ ■ ■ E ^.-.0'. *i,---, *»-o^ • • • «fc' 



• r exp (-/ £ X, +m (l - xj) dfl?'(0 

•'0 \ m=l / 

oo oo / /JV-1 

+ 2 ••■ Z(i-<2>, 



• fi(njv)Pf-i(0, n! , • • • , ny-Os?' • • • x"/ 
+ T E P*(0, • • • , 0) f exp (-t E X, + «(l - <)) dB*(0 

(t = 0,1, .-• ,N). (4) 

The integrals on the right side of (4) are recognized as the Laplace- 
Stieltjes transform (/?,•(•))' of the i-busy period distribution function 
with argument 2»-i A <+>»(1 — x m ). Hence (4) yields the set of simul- 
taneous functional equations 



9 



+ t ft(E Wi - *-)) E *vo, ■ ■ ' , 0) 

A \ m =i / 4-0 

-P.._ 1 (0 1 ••• ,0) (i = 0, 1, ••• ,N). (5) 
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V. SOLUTION OF THE FUNCTIONAL EQUATIONS 

For notational convenience, we define the nesting operator S for any 
sequence of functions {/*(■)} f° r which it is meaningful: 

*-0 

We shall denote by x the vector with components x x , • • • , x N , and 
by the vector with all components zero. Both vectors and vector- 
valued functions will be denoted by boldface type, and square brackets 
will be used to enclose vector arguments of vector-valued functions. 
Finally, we will denote by <£(v) the first component of a vector v. 
Define the vector functions 

Zi[x t , ■•• ,x N ] = I /3.|Z) A. +m (l - x m )J , x x , ■ ■ ■ , x N . x 

(i = 0,1, •■• ,N) (6) 
so that (5) can be rewritten 

<7,(x) = (A-.CZ.-M) + ^(Z,[x]) f;P*(0) - P,-,(0) 

A * = 

(i = 0,1, ... ,N). (7) 
Iterating v — 1 times on i in (7) we obtain 

gi (x) = g t J'z Z,_Jx]) + I fp t (0) £ X,_ m </ S Z,_ t [x]) 

Vt-0 / A £ = o m = \* = / 

- I)Pi-»(0) (i = 0,1, ■■■ ,N;p = 1,2, •••)• (8) 

m-l 

In particular, when v = N + 1 (8) can be written 

?l .(J; z,_ t [x]) - ^,(x) = p(o)(i - i g x -^(J o z -*w)) 

(t = 0,1, ••• ,N) (9) 
where we have set 



V 

I 

A-0 



P(0) = 1^(0). (10) 



We shall now solve (9) by extending a method devised by M. F. 
Neuts 12 for the solution of a related equation in one variable. 
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Define the iteration procedure 
V, ( ' } [x] = 8 Z l ._JV. ( '- 1) [x]] 

k-0 

(t-0,1, ••• ,N;Jr 1,2,... ;VrW = *)• (11) 
Using (11) in (9) gives 

<7.(v:"[x]) - *art ,_l) M) 

= p(o)(i - 7 i: xi.-«( s z,._ fc [v, o -"[x]])) 

\ A m = o \fc=o ' ' 

(1 = 0,1, •■• ,N;j = 1,2, ■••). (12) 
Adding equations (12) for ; = 1, 2, . . . , n yields 

</.(vHx]) - gM = P(o) £ (i - * £ x.-^h z.-*[Vi»M])) 

(t = 0, 1, ••- ,JV;n = 1,2, •••)• (13) 
Now let n -> oo in (13). We will show in the next section that 
lim VHx] - 1 (ft S! 1. • " • » *# i 1; » - 0, 1, ■ ■ • , JV) (14) 

n-«oo 

where 1 = [1,1, • • ■ , 1], so that (13) becomes 

*(i) - g<to = P(o) £(i-{£ a.-J s z.._jv^ ) [x]])) 

(i = 0,1, ••• ,N). (15) 
Notice that 

£ g&) = 1 (16) 

1-0 

and 

£ gM = P(0) (17) 

so that upon setting x = in (15) and adding for i = 0, 1, ■ • • , N 

we obtain 

p(o) = (i + £ Am) ' (is) 
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where 
A,(x) = E (l - I Z h-J 8 Z,_ t [V,"W 

(t = 0,1, •■• ,N). (19) 

(We remark that P(0) ^ 1 - 2f-o M* because the set of switch points 
is not an arbitrary subset of the set of all points at which units leave 
the server.) 

It remains to calculate </,(l). Physically, 0,(1) is the probability that 
at the instant the server leaves some queue, that queue is queue i. This 
event occurs if 

(0 the last time the server left a queue the system was empty, and 

the next arrival occurred at queue i, or 
(it) the last time the server left a queue the system was not empty, 
and the queue was queue i — 1. 
Event (i) has probability (X,//\)P(0); event (ii) has probability £7.-i(l) - 
gr<_i(0). Hence 

0,(1) = ^P(O) + (^.(l) - f7.-i(0)) (t = 0, 1, -•• , N). (20) 

[Equation (20) can also be obtained directly from (7) with x = 1.] 
But the difference (ft-i(l) - <7,-i(0)) can be evaluated from (15) with 
x = 0. Hence 

0,(1) = h.p( ) + P(0)/1,._,(0) (i = 0,1, ... , N) (21) 

so that (15) can be rewritten 

^ + ^,._,(0) - A,(x) 
g,(x) = ± - N (i - 0, 1, • • ■ , N). (22) 

1 + £ AM 

The quantities on the right side of (22) are completely specified; the 
set of simultaneous functional equations (5) has been solved in the 
sense that g t (x) may be calculated for any x fS 1. 

VI. PROOF OF CONVERGENCE 

We wish to prove statement (14) : 
lim V{ B) W = 1 (z, £ 1, • • • , x H £ 1; i = 0, 1, • • • , N). 
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Note first that V| B> [x] is a vector whose (N + 1 - ro)th element 
(m = 1, 2, • • • , N) is 

is:z,._ i [vr ,) W])- 

Therefore, we need show only that 

timJ* z.-avj-'w]) = 1 



(i = 0, 1, • ■ • ,N;m= 1,2, ■•■ , N; Xi £ 1, • • ■ , x„ 5 1). (23) 

From the definition (11) it is clear that the sequence { V\ n) [x] ) is bounded 
as n -> oo for x ^ 1, and therefore the sequence {<7,(V' n> [x]) } is bounded 
as n — * oo for x ^ 1. Also, 

(24) 



k-0 



0^ E Z.-JVrix]] ^ 1 (n > 1, x ^ 1). 



We now turn our attention to equation (13). From (24) we see that the 
right side of (13) increases monotonically with n for x ^ 1, and there- 
fore the sequence (ff,(V- n) [x]) } increases monotonically with n for x ^ 1. 
Thus the sequence |fl'.-(Vi" ) [x]) } is monotonically increasing and bounded 
for x ^ 1, and therefore has a limit. Hence the left side of (13) has a 
limit, which implies that the series of nonnegative terms on the right 
side of (13) converges. This in turn implies that 

lim \ f) \,-„*( H Z,_ t [V! n, [x]]) = 1 (x =g 1). (25) 

n— oo A m-0 \* = ' 

Statements (24) and (25) together imply (23) , completing the proof. 

VII. MEAN NUMBERS OF WAITING UNITS 

Denote by n.-(fc) the mean number of units waiting in queue i + k when 
the server leaves queue i(i = 0,1, ••■ ,N;k = 0,1, ■•• ,N; n,(0) = 0). 
For convenience, let fr(l)fl<(fc) = fn t (k) and m.-(l) = m, . Then £?-o ™< 
is the mean number of waiting units found by the server in the next 
queue in cyclic order at a switch point, and 2Jt~o ™<( fc ~~ ■) * s *^ e 
mean number of waiting units in queue k at a switch point. We shall 
evaluate m,(/c) (i = 0, 1, ■ • ■ , N; k = 1, 2, ■ ■ ■ , N). 

We first note that m<(fc) is given by 



ox* 



' , Xir) 

(i = 0, 1 



,#;* = 1,2, ••• ,N) (26) 
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and the mean 1-busy period b. ( = ft</(l — AA) generated by a unit 
in queue i is given by 



*--!*« 



(i = 0, 1, ■•■ ,2V). (27) 

Differentiating through (5) we obtain 



£ *<* .••■.**) 



0<(£) X, + .(l - xS) jjf fl<-i(fc , *i , • ' ' , 



+ (1 - «(JV— *)) ~ p,_,09, ,.r,,---, a*-i) 

(i = 0,1, •■■ ,N;k = 1,2, ■■■ ,N) (28) 

which upon setting X\ = . . . = .r.v = 1 gives the two-dimensional set 
of linear equations 

m^k) = X, +t 6,«,-., + ^P(0)X I+ *&, + (1 - «(#-*)) >»<-i(* + 1) 

A 

(i = 0, 1, ■ ■ • , N; k = 1, 2, • • • , N). (29) 

For each i, (29) can be solved successively starting with k — N 
and working backward : 

y + 1 i+i 

ro,(AT - j) = X, + iV _, £ 6.- + i-,w,-_„ + X _1 P(0)X, +iV -y Z X, +1 _,&, + .-, 
►=i »=i 

(« = 0, 1, ••• ,AT;j = 0, 1, ••• ,N - 1). (30) 

In particular, when j = N — 1 equation (30) can be written 

i+JV i + .V 

M t = X 1 + 1 £ & r + 1 ?n. r + X-'P(0)X.--i E Xr + i6r + i 

v - i i + 1 r - 1 + 1 

(t = 0, 1, ■•• ,JV). (31) 

When X, + i6 1+l m, is added to both sides of the ith equation of the set 
(31) we have after rearrangement 

i+JV i + JV 

01,(1 + X, + 1 b, tl )X7'. - X-'P(O) £ K + ib, + i = E b, +1 m, 

►=i+l y-i 

(i = 0, 1, ••• , N). (32) 
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The sum on the right side of (32) is a constant independent of the 
value of the index i. Hence 



?.,(1 + X 1+1 b, + >)X~ + \ - X" l P(0) £ X„ +1 6, +I 

v = i + 1 

= m, + ,(l + X,- +J+1 &, +I+I )X, 7 + 1 / + i - X-P(O) £ X, +I 6, + , 



p = i + 1 

i + j + JV 

E 

r-«+»+l 

(i = 0, 1, ••• ,iV;i = 0, 1, ••■ ,N). (33) 
Combining (33) and (31) yields 



m . = ^tip( ) p Pi+1 (t = 0, 1, • • • , N) (34) 

A 1 — p 

where we define p, = X t 7i,- and p = XXo Pi • Note that for (34) to be 
meaningful we must have p < 1. The {m,-(fc)l can now be calculated 
from equations (34) and (30). 

VIII. LAPLACE-STIELTJES TRANSFORM OF CYCLE TIME DISTRIBUTION 
FUNCTION 

Consider the set of switch points associated with queue i, and append 
to this set every switch point associated with queue i — 1 at which the 
server finds the system completely empty. Call the elements of this 
augmented set the record points associated with queue i. 

We define the partial cycle time for queue i as the elapsed time be- 
tween a switch point associated with queue i — 1 and the temporally 
preceding record point associated with queue i. Denote by (?,•(•) the 
distribution function of the partial cycle time for the i'th queue, and by 
■?,(•) its Laplace-Stieltjes transform. 

Since queue i is necessarily empty at an associated record point, all of 
the units waiting for service in queue i at a switch point of queue i — 1 
must have arrived during the preceding partial cycle time. Let P,-i0") 
be the conditional probability that j ^ units will be waiting for 
service in queue i, given that a switch point associated with queue i — 1 
has just occurred. Then the distribution function (?<(•) of the partial 
cycle time for the ith queue and the distribution {P t -_i0")} of the 
number of units that arrive (according to the Poisson process with rate 
X.) during the partial cycle time are related as follows: 

P.-.0') = f^# exp (-X.-0 <*&(Q 

Jo J ■ 

(i = 0,1, ••• ,N;j = 0,1. •••)• (35) 
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Notice also that the distribution {P,_i(j)} has probability generat- 
ing function 

E P«-i&V = 0i - l(X ' h ,\,' — (i - 0. 1. ■ • ■ , JV). (36) 

Substitution of (35) into (36) yields, for the Laplace-Stieltjes transform 
-/,(•) of the partial cycle time distribution function for queue i, 

Us) = K -^ 77T (i - 0, 1, •• • , JV). (37) 

We define the (full) cycle time for queue i as the partial cycle time 
plus the time required to serve those units, if any, waiting in queue 
i when the server finishes queue i — 1. (Notice that in order to be 
counted as a cycle for queue i, a time interval must contain a partial 
cycle ending at a switch point at queue i — 1.) Denote by (?,•(•) the 
distribution function of the cycle time for the ith queue, and by *y«(0 
its Laplace-Stieltjes transform. 

The cycle time distribution function (?,(•) is related to the partial 
cycle time distribution function (?,(•) as follows: 

G«(0 = f X. %7- ex P (-U)Bf(t - & d6 t ® 

■' O f - J ■ 

(Bf(-) = l;t = 0,1, ■■■ ,A0- (38) 
Taking Laplace-Stieltjes transforms throughout (38) we obtain 

Y.-(s) = 7.(X.- + s - X,-/S.(s)) (» = 0, 1, • • • , N). (39) 

Hence we have for the Laplace-Stieltjes transform yi(s) of the cycle 
time distribution function for the ith queue 

( MM - g i ^ 

7,(S) = - — 77T d = 0, 1, • • • , JV). (40) 

By differentiating through (40) we obtain for the mean cycle time 
i; the intuitively obvious result 

I. = (&. + X" 1 )^,.., (t = 0, 1, • • • ,N). (41) 

IX. THE GATING MODEL 

Consider now a system of N ^ 1 cyclic queues described by the 
gating model of Section I. Define P,(n, , • ■ • , n N ) as the joint prob- 
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ability that at the instant the server leaves a queue, that queue is 
queue i (i = 0, 1, ■ ■ ■ , N — 1) and th units are waiting in queue i + 1, 
?i 2 units iD queue * + 2, • • • , and n N units are waiting in queue i (that 
is, n N units arrived at queue i after the closing of the gate). Denote by 
#*'(•) the ;-fold convolution with itself of the service time distribution 
function #,(•)• Then 
P,(n, , • • • ,n N ) 

= 2Z- Z P«-i(i, *i . ■ ■ ■ , *»-0 f ri &+«(»« - *. ; ') 

•Q, +JV (n„ ; /) dHf(i) + P,--,(0, n, , • • • , n w -,)(l - s( Z ». 



• «(%) + T Z Pk(0, ■■■ ,0) T fl Q i+m (n m ; dHM 

A i=o J m=l 

(i = 0, 1, ••• ,N - 1) (42) 

where Q t = Qi +A r. 

Equation (42) , for the JV-queue gating model, is only trivially dif- 
ferent from (1), which describes the (N + 1) — queue exhaustive ser- 
vice model. The analogue of (5) is 



g { (xi , ■ • ■ , x N ) = 0«-iU<( Z) X I+m (l - x m )j , x t , ••• , x K -i 

\ ( N \N-l 

+ Y H Z Wl - «-) ) Z P*(0. ■ • • , 0) 

-P,-,(0, ■•• ,0) (t = 0,l, ■•• ,iV-l) (43) 

where ^{(O is the Laplace-Sticltjes transform of the distribution func- 
tion Hi(-), and (/{(rr-i, . . . , x N ) is now the generating function for the 
gating model state probabilities. The solution of (43) follows that 
given for (5), and a complete analysis may now be carried out in a 
manner similar to that employed for the exhaustive service model. 
(We remark in passing that the equations originally considered by 
Neuts are those of the gating model with N = 1.) 

X. SUMMARY 

Two models of a system of queues served in cyclic order by a single 
server have been presented. One of these, the exhaustive service model, 
has been analyzed in detail. This model is described by the imbedded 
Markov chain probability state equations (1), from which a set 
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of functional equations (5) for the probability generating functions 
are derived. The functional equations are solved with the help of a 
generalization of an iteration procedure used by Neuts. The equa- 
tions (5) are then used to obtain explicit expressions for various 
mean values, such as the mean number of units found waiting by the 
server in the ith queue, given by equation (34), and the mean cycle 
time, given by equation (41). The Laplace-Stieltjes transform of the 
cycle time distribution function is given, in a form suitable for nu- 
merical computation (and hence numerical inversion), by equation 
(40). 

It is then shown that the gating model is described by state equations 
only trivially different from those of the exhaustive service model. It 
is now easy to adapt the methods and results of the detailed analysis of 
the exhaustive service model to a similar analysis of the gating model. 

It is noteworthy that all results are expressed directly in terms of the 
single state probability P(0) and the relevant generating functions, so 
that there is no need to evaluate the individual state probabilities. The 
calculations are thus reduced to the iteration algorithm, which may be 
suited to digital computer solution. 
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